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In this work, one shows that a specific non-minimal coupling between the scalar curvature and 
matter can mimic the dark matter component of relaxed galaxy clusters. For this purpose, one 
assesses the Abell Cluster A586, a massive strong-lensing nearby relaxed cluster of galaxies in virial 
equilibrium, where direct mass estimates are possible. The total density, which generally follows a 
cusped profile and reveals a very small baryonic component, can be effectively described within this 
framework. 
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I. INTRODUCTION 

In the last two decades, cosmology has acquired a new 
standing due to a wide range of observations of high pre- 
cision data [3 12] . These observations seem to reveal that 
the two dominant components of the universe are a non- 
baryonic form of matter, the so-called dark matter [3], 
and an exotic form of energy, dark energy [4] . The former 
leads to structure formation and the last is responsible for 
the present accelerated expansion of the universe. The 
search for an explanation for the existence and properties 
of these dark components of the Universe has prompted 
a strong debate about their origin and nature. 

The behaviour of the galaxy rotation curves and the 
mass difference in clusters of galaxies suggest the exis- 
tence of a non-interacting form of matter, the so-called 
dark matter, at galactic and extra-galactic scales respec- 
tively. The presence of this sort of matter is revealed 
only through its gravitational interaction, since it pre- 
sumably does not interact via any other detectable way. 
For instance, an evidence for the existence of dark matter 
is the velocity dispersion of clusters that, combined with 
its morphological features, leads one to conclude that the 
overall mass should be far greater than the visible mass. 

Amongst modified gravity theories, those which intro- 
duce higher order curvature terms in the Einstein- Hilbert 
action — and in particular the f{R) theories of gravity 
[51 [B] — offer alternative explanations for these problems, 
along with other strong motivations coming from cosmol- 
ogy, astrophysics and high-energy physics 7 . 
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Recently it was shown that, in the context of power- 
law theories, f{R) ^ i?", the rotation curves can be ex- 
plained as a curvature effect [SI [9] ; a similar result [10] is 
also obtained in another extension of gravity that relies 
not only on a non-trivial curvature term, but also on a 
non- minimal coupling between matter and geometry [llj . 

Clusters are the largest structures whose masses can 
be measured in a reliable form, since they constitute the 
largest configuration of objects that passed through grav- 
itational relaxation and entered into virial equilibrium 
[12) . The measuring methods. X-rays and gravitational 
lensing for visible and dark matter respectively, reveal 
that the masses of clusters are approximately seven to ten 
times larger then the total combined mass of stars and 
hot gas in the intracluster medium (ICM). This missing 
mass is predominant in the inner regions of clusters, but 
it is also extended beyond the core radius, as defined by 
the gas density distribution. 

In the context of modified Newtonian Dynamics 
(MOND) [13], another popular (albeit incomplete [TJ]) 
alternative to General Relativity (GR), it is difficult to 
explain cluster observations without the assumption of 
"real" dark matter, because the mass of the dark compo- 
nent is otherwise significantly reduced but not completely 
removed [15) . 

The problem was recently introduced in the context 
of f{R) theories, by considering a generalized version of 
the virial theorem [16) . A similar proposal, based upon 
a possible unification and interaction of dark matter and 
dark energy, had been put forward before, using similar 
considerations [XT'. Assuming a steady state, the virial 
theorem can be used to deduce the mean density of as- 
trophysical objects, and therefore the corresponding to- 
tal mass, such as galaxies and clusters, by observing the 
velocities of test particles rotating around them. The 
overall mass of a sample of clusters can also be estimated 
through corrections in the gravitational potential that 
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emerge in the weak field limit of these theories . Along 
with these works, there were also conducted the first clus- 
ter abundance constraints on a modified gravity model, 
specifically the modified action f{R) model, which pre- 
dict deviations in the abundance of massive dark matter 
halos HH. 

Following Ref . (TTj , it has been argued in several works 
that both dark energy and dark matter can be effectively 
described within the framework of a non-minimal cou- 
pling between geometry and matter, together with other 
cosmological and astrophysical implications: these in- 
clude the perturbed hydrostatic equilibrium of the sun 
[5D] , the acceleration of the universe |[2J| , the galaxy ro- 
tation curves p3], the reheating scenario after infiation 
, the mimicking of a Cosmological Constant [13] and 
the change in the gravitational potential 

In order to extend this non-minimally coupled model 
to the study of clusters, one explores here the scenario 
of the Abell cluster A586. This cluster is a nearby mas- 
sive strong lensing cluster, with a high X-ray luminos- 
ity, related with the mass of the intracluster gas. Its 
morphology is well suited to the assumption of spherical 
symmetry, making it an ideal testbed for modifications 
of gravity at extra-galactic scales. 

This work is organized as follows. First, one shortly 
reviews the formalism leading to the total cluster mass 
and the field equations that result with the presence of 
the non-minimal coupling, and sees how to work out the 
problem of mimicking the dark matter component of a 
cluster. Then one studies the particular case of the Abell 
cluster A586 and discusses the coupling function that 
leads to the dark matter mimicking. In the last sections, 
the numerical results are presented and discussed. 



II. GALAXY CLUSTERS 
A. Hydrostatic equilibrium equation 

Galaxy clusters are usually considered as closed gravi- 
tational systems with spherical symmetry and in hydro- 
static equilibrium when virial equilibrium is attained — 
despite the fact that recent observations reveal that clus- 
ters in general have more evolved structures with strong 
interactions and dynamical activity, particularly in their 
inner regions [55', "55] . 

Under these assumptions, the structure equation can 
be derived from the coUisionless Boltzmann equation for 
an isotropic system [18l [27] 



^ {Pg <^l) - -Pg 



dr 



(1) 



where $ is the gravitational potential of the cluster, (Jr 
the mass- weighted velocity dispersion in the radial direc- 
tion, and Pg the gas mass-density. The pressure profile 
is related with these quantities through, P = cr'^Pg, and 
inserting in Eq. ([l]) leads to P' = —pg^', where the 



prime indicates differentiation with respect to the radial 
coordinate. 

For a gas sphere, with temperature profile T(r), the 
velocity dispersion becomes ct^ = J^i < > —vf — 
kT/pTTip, where k is the Boltzmann constant, p « 0.609 
the mean mass particle, and mp the proton mass. Intro- 
ducing this expression into Eq. ([l]) gives 
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Therefore, from the models for the gas density and 
temperature profiles, the total mass enclosed in a sphere 
with given radius r can be obtained from the hydrostatic 
equilibrium equation 
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and is composed by different contributions: gas mass -I- 
mass of galaxies -I- cD-galaxy mass + dark matter. In 
general, the galaxy part contribution is considered negli- 
gible with respect to the other two components such that 
M ^Mg + Mdm- 

The equilibrium equation, Eq. Q, which leads to the 
total mass within a cluster, can be used to derive the 
amount of dark matter present in a cluster of galaxies, 
and its spatial distribution, by the mass difference be- 
tween this and the gas-mass estimates that are provided 
by X-ray observations. 

In the following, one tries to reproduce the total mass 
profiles for a particular cluster, the Abell cluster A586, 
where the additional density component emerges not 
from an exotic form of "dark" matter, but from a change 
in geometry due to the presence of the non-minimal cou- 
pling in the action. 



B. Cluster density profiles 

Clusters of galaxies present a velocity dispersion with 
a relative constant value along the radius, which reveals 
an implicit matter density profile of the form pM cx: r~^. 
The singular isothermal sphere is the simplest model 
whose density profile, pM = (Jy/2nGr'^, leads to a con- 
stant and isotropic velocity dispersion along the clus- 
ters profile. 

However, numerical simulations to the observational 
data show that the density profile of clusters are flatter 
than the isothermal one at small radii and steeper at 
larger radii. A general form replicating this behaviour is 
given by the cusped density profile. 
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Pm = Po 



(5) 



where the parameters p and q determine the respective 
slopes of the power-law density profile in the inner and 
outer regions, and is the so-called core radius, signaling 
the steepening of the profile. These parameters usually 
present values within the ranges 1 ^ p ^ 1.5 and 2.5 ^ 
p < 3, for models such as the NFW profile, for p = 1 and 
q — 3, the Moore profile, for p = 1.5 and g = 3, and the 
Rasia profile, for p = 1 and q = 2.5. In fact, one of the 
important results of hierarchical CDM models is that the 
density profile of cold dark matter halos is described by 
the Navarro- Frenk- White (NFW) model p8] . 

The more relevant baryonic component is the ICM in 
the form of hot gas, which emits observable X-ray radi- 
ation. In general, the observed surface-brightness pro- 
files of clusters, that reveal the underlying visible matter 
profile of this intracluster gas, are well described by the 

so-called /3-model of the form pg = po (l + f'^ 1^1) '^^^^i 
where Tc is the core radius, which has a specific value for 
each cluster [29l. 



C. Abell Cluster A586 

In this cluster, the emission measured profile is de- 
scribed by the following model 



Pg{r) = P90 1 + 



,2 \ -3/J/2 



(6) 



From a least-square fit of the available data, the obtained 
values for the free parameters are (3 = 0.518 ± 0.006 and 

= 67 ± 2 h^Q kpc [3D| . This analytical expression re- 
sults from a simple modification of the usual /3-model. 
The ICM particle number density profile is given directly 
by the analytic fit to the projected emission measure pro- 
file, Eq. ([6]), and is easily converted into the gas density. 

The deep potential well of the cluster compresses the 
ICM gas to X-ray emitting temperatures, such that the 
gas temperature can be inferred from the derived spec- 
trum. For a polytropic temperature model [50] . 



T{r) = To 1 + 



-3/3(7-l)/2 



(7) 



where and /? are the parameters obtained in the fit 
for the visible mass density, a least-square fit leads to 
the values Tq = 8.99 ± 0.34 keV and 7 = 1.10 ± 0.03 
|30j . Naturally, setting 7 = 1 leads to the isothermal 
temperature profile. 

The total mass follows from the hydrostatic equilib- 
rium equation, Eq. Q, using the previous expressions 
for the visible density and temperature profile 



M{r) = Mo 



(8) 



with 



Mn 



Gprup 



(9) 



One may approximate the dark matter density by 
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where 



PdmO 



Mo 



(11) 



and one defines e = [1 + 3(1 — 7)]/?, for convenience. 

Finally, one introduces the virial radius ry, which sig- 
nals the distance where the density of the cluster ap- 
proaches a multiple a of the critical density of the uni- 
verse pc = 3Ho/8ttG, where Hq is the Hubble parameter 
in the present time. Thus, from pdmirv) = ctPc one may 
estimate the virial radius from the relation Eq. (10 1, 



PdmO e 

Pc a 



(12) 



which, for pdm ~ 200pc, with a = 200, leads to a virial 
radius of ry ~ 1.497 Mpc. 



III. THE NON-MINIMAL COUPLING 

A. Field equations 

Recently, it was considered a generalization of the f{R) 
modified theories of gravity that further extends the pres- 
ence of curvature invariant terms in the Einstein- Hilbert 
action by coupling non-minimally the scalar curvature 
with matter , 1 Ij . The action is obtained by adding an 
additional f{R) term, dependent on the scalar curvature 
R and coupled with the matter Lagrangian Cm, to the 
usual metric theories of gravity, such that [11] 



S 



-h{R) + [l + f2iR)]C„ 



~gd X 



(13) 



where fi{R) (with i = 1, 2) are arbitrary functions of the 
scalar curvature and g is the metric determinant. Setting 
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f2{R) = one obtains the usual f{R) theories of grav- 
ity, and GR is recovered if one takes the hnear function 
/i(i?) = 2kR, where k = c^/IGttG. 

Variation with respect to the metric, g^jy, leads to the 
field equations. 



Uf^ the four- velocity (with Ufilf^ = —1), the respective 
trace is, T = 3p — p, such that one may write. 



2£ = 3p~ p 



p^{2C + p)/S 



(19) 



{Fi + 2i^2^m) Rfj.u — 2 /i5m<^ — 

=- (Dmi^ - 9,..n) {Fi + 2F2£m) + (1 + /2) Tf,. , (14) 

where is the energy-momentum tensor. The respec- 
tive trace equation is given by 



(Fi + 2F2£„0 R-2h = 

= -in [Fi + 2F2C^) + (1 + h) T 



(15) 



B. Non-conservation of the energy momentum 
tensor 



However, one obtains different equations of state de- 
pending on the chosen Lagrangian density £. For in- 
stance, ii C ~ —p |31j . one has the state equation 
p = —p/3, a perfect fluid with negative pressure, not 
a cosmological constant with p\ ~ —ph_\ and for C = p, 
then p — p, which is the equation of state for ultra-stiff 
matter. Therefore, this non-conservation is in fact a fun- 



damental property of the model, Eq. ( 13 1, meaning that 



even under a suitable conformal transformation the en- 
ergy momentum tensor is not conserved for all different 
types of matter. 

Assuming that matter is described by a perfect fluid, 
the non-conservation, Eq. (16), results in an extra force 
of the form 



From the covariant derivative of the field equations, 
Eqs. (14), along with the Bianchi identities, V^G^j^ = 
0, one encounters the covariant non-conservation of the 
energy-momentum tensor |11| 



T^,)V^log(l 



which can be regarded as an energy-momentum exchange 
between matter and geometry. The usual covariant 
conservation for the energy-momentum tensor is natu- 
rally obtained when the non-minimal coupling vanishes, 

/2(i?)=0. 

This non-conservation implies that a massive particle 
in the absence of forces will not describe a geodesic curve. 
Strong variations on the extra term in the last equation 
can imply a violation of the Equivalence Principle, lead- 
ing to a possible way of testing and setting bounds on 
the coupling functions, f2{R)- Contrary to the usual 
Jordan-Brans-Dicke theories, it is impossible to perform 
a conformal transformation to the Einstein frame such 
that the coupling disappears and the covariant conserva- 
tion of the energy-momentum tensor is recovered for all 
matter forms [35]. 

Indeed (following Ref. (22), using the conformal trans- 
formation 



p+p 



{C^ + p) log (1 + J2) + VuP 



hf"'{20) 



which has units of an acceleration, and where the projec- 
tion operator is given by ft,^^ = g^m — U^Uv^ such that the 
J2) , (16) extra force is orthogonal to the four-velocity, h^^U^^ = 0. 



C. Mimicking the dark matter component 

In order to test the effect of the non-minimal coupling, 
one focuses on possible deviations from GR due to the 
presence of the coupling by considering the linear func- 
tion fi{R) = 2kR; following Ref. [10], one adopts a 
power-law non-minimal coupling 



/2 (R) = 



R 

Rn 



(21) 



where i?„ is a characteristic curvature. 

Considering a pressureless perfect fluid, Tfj,^ — pUfj,U^, 
and the respective trace equation, T = — p, along with 
the Lagrangian density, Cm = ~P choice discussed in 
depth in Ref. [3T]), the fleld Eqs. (Il4| become 



ffpi. -> Sa"!^ = h9p.i^ : (17) 
one obtains the covariant conservation law 



(18) 



only if Ti"' = f-^Ti''' and 2C = T. From the energy- 
momentum tensor of a perfect fluid, T^^ = {p+p)u^u^ + 
pg^u, where p is the energy density, p the pressure and 



n f R 
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1 

"2k 
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Rn 



and the respective trace, Eq. (15), reads 



(22) 
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R = 



2k 



1 + (1 - 2n) 



R 

Rn 



3n 

p □ 

K 



R 

Rn 



(23) 

If one considers a strong coupling (R/Rn)" » 1, the 
above has an imphcit, "static", solution, given by the 
vanishing of the second term in the r.h.s., 



R 



Pdm 
Pn 



(1 - 2n)- 



Pn 



1/(1- 



(24) 



where one introduces the characteristic density, p„ = 
Considering the tensorial character of the field Eqs. 



(14), one obtains the following relation between the en- 



ergy density and pressure of the mimicked dark matter 
component [TU] . 



Pdn 



Pdm = 



1 - n 

I -An 
n 

1 - An' 



2nR , 
2kR . 



(25) 
(26) 



From Eq. (24), one has a relation between the visible 



and dark matter density profiles, 



Pdn 
Pn 



1-n 
1 - 4n 



(1 - 2n) 



1/(1-") 



(27) 



The above clearly shows that one may obtain a "dark" 
component from the visible matter density profile, with a 
change of slope given by the scaling exponent 1/(1 — n). 
For convenience, Eq (23), can be rewritten in a dimen- 



sionless form, by defining the rescaled quantities |TU], 



D. Non-geodesic motion 

Before dwelling into the analytical and numerical de- 
tails of the scenario under scrutiny, one must first con- 
sider a possible issue: the presence of an extra-force, as 
given by Eq. (20 1. Considering a pressureless perfect 
fiuid along with the Lagrangian, Cm = —p, the same 
equation reduces to 

/^ = -V,log(l + /2)/i^'' . (32) 
For spherical symmetry, the radial component is 



r = -9,iog(i + /2) 



F2{R) 
l + f2{R) 



R! 



(33) 



where the prime denotes derivative with respect to the 
radial coordinate. Considering a power-law coupling of 
the form of Eq. (21 ), one has 



r 



(R/Rn)" R' 
''l + {R/RnTli 



(34) 



If the curvature R decreases along the radius, a negative 
exponent n will produce an inward binding force, and 
a positive one an outward force, which could break the 
stability of an orbit. 

As already seen in a previous work [TUj, different so- 
lutions for the differential equation, Eq. (31), lead to 



different behaviours of this extra-force, depending on 
the relevance of the coupling function. If it dominates, 
(i?/i?„)" » 1, then the force becomes 



r 

P_ 

Pn 

2kR 

Pnd 



1/(1-"') 

1 R 

9 Rn 



where one defines the length scale, r, 
resulting dimensionless form of the trace Eq. ( 23 ) is 



(28) 
(29) 
(30) 

l/^fR^,. The 



+ (1 - 2n)^^" 



(31) 



where □ 



y 



'^(P/d?y. The equivalent to the "static" 



solution, Eq. (24), is thus, g = {I — 2n)^/^^ "\ as may 



be easily checked. 

Concluding, one uses the visible matter density profile 
p, which in a cluster is approximately equal to the ICM 
gas density, p fa pg, in Eq. ( 29 ) , to solve the differential 



equation Eq. (31), leading to a dark matter component 



that can account for the large mass difference in the clus- 
ter total density profile, as will be discussed in section 

ttvH 



fr ~ -n 



R 



(35) 



On the contrary, a perturbative regime leads to a sup- 
pression of the expression above by a factor (i?/i?„)" <C 
1, 



R 

Rn 



R 



(36) 



The "static" solution given by Eq. ( p4| is derived from 
the dominant condition f2{R) 3> 1; taking only the outer 
slope behaviour of cusped density profiles for the dark 
matter component, pdm oc r"™, and a positive exponent 
n (as will be used in the following section), one concludes 
that the extra force has an outward direction. 



R' 

'r 



Pd7n 



nm 
r 



(37) 



For n and m of order unity, one obtains a force 
much greater than the Newtonian counterpart, /^r = 



6 



GM{r) /r'^ , thus destabilizing the orbital motion of test 
particles. In the inner region, r <C Tc, one has an ap- 
proximately constant dark matter density, pdm ~ PdmOi 
yielding a negligible extra force. 

However, it was shown in Ref. [10, that the most 
general solution of the differential, Eq. (31), is not the 
"static" solution, g — {1 — 2n)^/(^~"\ because the gradi- 
ent terms on the r.h.s. of this equation are dominant: the 
system instead follows a "dynamical" solution, g, that ex- 
hibits very small oscillations around a constant value. It 
was shown that, although the above analytical treatment 
remains valid, this implicates that the non-minimal cou- 
pling is actually perturbative, i.e. /2 — (R/Rn)" ^ 1- 

For a fixed it is then clear that f2{R) increases for 
lower values of a positive exponent n; thus, one expects 
that there is a lower bound on n, such that the condition 
< fif fails due to the insufficient suppression of the 
extra force /''. 

In order to assess the value of this lower bound, one 
first uses Eqs. ^ and ( 10 ) to write the following approx- 
imations, valid in the large radius region r ^ re'. 



Mir) - Mo 



Pdm ~ PdmQf- 



(38) 



-2+e-/3 



Using the above, one may write 



Vi?'V R \ Pn J Pdm 



(39) 



(-2+e-/3)Ti-l 



so that, writing pdmo — 2K/rQ, the perturbative condition 
<^fj^ = GM{r)/r^ reads 



Thus, a lower exponent n will increase the l.h.s., as 
expected. The following section will discuss this lower 
bound on n in the case of the Abell cluster A586. 



IV. THE ABELL CLUSTER A586 CASE 
A. Density Profiles 

After considering the possibilities offered by the non- 
minimally coupled model in the previous sections, one 
now implements the mimicking of the dark matter com- 
ponent of the cluster in consideration. 

Recalling Eqs. ^ and (38), one writes the approxi- 
mations valid in the large radius regime r ^ Tc, 



Pa 



-3/3 



: Pdm 



-2-3(7-l)/3 



(43) 



Therefore, from the scaling given in Eq. ( 27 ) , one obtains 
the relation pd„i ^ P^^^ allowing for an estimation of 
the expected value of the n exponent: 



2 + 3(7-2)/3 

2 + 3(7- l)/3 "-'^^ ' 



(44) 



using f3 = 0.518 and 7 = 1.10 30 ; and, for the isothermal 
case with 7 = 1, 



1 



3/3 



(2 - 7) w 0.223 



(45) 



Hence, the results coming from the presence of the non- 
minimal coupling should be similar of whether one is 
considering an isothermal or a polytropic temperature 
profile. 



v/2n(2-e + /3)^ f ^J" ^ ' ^ 



ro / \r. 



< 1 



(40) 

One now assumes that — 2n + (e — /3)(n — 1) > (as 
is the case with the Abell cluster A586, discussed in the 
following section). Thus, the relevance of the extra force 
decreases with the distance from the cluster's center, and 
one may saturate the l.h.s. of the above equation by 
setting r ~ Tc, yielding 



6(r„,n) = v/M2^7T^fVe-) -«1 ■ (41) 

Finally, notice that, for a positive exponent n, the per- 
turbative condition /2 <C 1 translates into 



Rn / 



< 1 ^ PdmO < 2Ki?„ ^ ro > r„ . (42) 



B. Numerical results 

In order to mimic the dark matter component of this 
cluster, one numerically solves the differential Eq. (31), 
for an assumed gas density profile, as given by Eq. ^), 
and varies the parameters i?„ and n in order to obtain 
a best fit between the derived mimicked dark matter 
density ptotai — Pgas and the observed curve. The con- 
sidered range for the numerical integration lies between 
ri = lOr^, the Scharwzschild radius of the Cluster, and 
a sufficiently large radius, above the observed size of the 
cluster. 

The results of the numerical simulations are shown in 
Fig. (II]) and the resulting integrated mass radial pro- 
file, Mir) = p{r')dr', in Fig. (B, for the values of 
n = 0.23, 0.33, 0.43, 0.53, and 0.637 with a best fit at- 
tained for n = 0.43; r„ = 10~^ pc is kept constant, as its 
value does not change the mimicked dark matter density 
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FIG. 1. Mass density profiles for the visible matter component 
(gray), the total density profile fi-om the hydrostatic equi- 
librium equation (dark), and the mimicked density profiles 
(dashed) for different values of n, where the dashing length 
increases with n. The best-fit obtained is for n — 0.43 (dark 
dashed) and the horizontal dotted line is for p = 200pc. 



TABLE I. Virial radius and respective mass for pdm ~ 200pc 
Mass differences for different values of the n evaluated at r = 
1.165 Mpc. 



profile (as discussed below). This value for the exponent 
n is determined by the minimization of the relative dif- 
ference between the obtained mimicked total mass and 
that expected from the hydrostatic estimative, Eqs. ([4| 
and (|8]), with values shown in Table 

As one can see, these results show that, in the outer 
regions of the cluster, the density decreases with the value 
of the n along with a decrease in the total mass for the 
cluster. The best fit n = 0.43, is the one that yields 
a mimicked profile which minimizes the mass difference 
to only 6.94% at ry = 1.165 Mpc, the virial radius for 
the density profile derived from Eqs. ^ and (|8|. In 
fact, in Fig. one can see that the values in the range, 
n < 0.43, lead to a dynamical mass larger than the upper 
limit set by the model, M = 4.25 x lO^Mg. 

Notice that this best fit n = 0.43 deviates from the 
expected value n 
curves for n 



0.279, computed in Eq. (|44|. As the 
0.33 and n = 0.23 show, this stems from 
the fact that the mimicked dark matter density arising 
from the expected value is unable to follow the observed 
profile. Indeed, while n = 0.279 yields the same slope 
as the observed total mass density, it leads to a right- 
shifted curve, giving rise to a much higher total mass for 
the cluster. 



FIG. 2. Mass radial profile, resulting from the integration of 
the mass density profile. The dotted line is for the dynamical 
mass of M = 4.25 X 10^" Mq (polytropic case) [30] . 
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FIG. 3. Resulting power-law couplings, f2{R) = (R/Rn)"', 
after the integration of the differential equation, Eq. (31 1, for 
different values of n. 



Numerically, one finds that the mimicked dark matter 
profile is almost independent of the characteristic length 
scale r„, provided one sets the bound < 1 Mpc. This 
independence is easily explained: a different r„ induces 



a change in the boundary conditions of Eq. (31), which 
translates into a change of the dimensionless solution g 
(not shown); however, this is counteracted by the pres- 
ence of r„ in the relation between the latter and the mim- 
icked dark matter density, Eq. ( 29 ) . 

Fig. ([4]) shows that values above r„ ~ 1 Mpc lead to 
a mimicked dark matter density that strongly deviates 
from the observed profile, with a sharp rise (and a sub- 
sequent fall into the expected outer slope). The onset of 
this rise decreases with higher values of r„ — i.e., lower 
values of i?„ . 

Despite the complexity of the differential Eq 



(311, 



the existence of an upper bound for this independence 
seems to be due to the fact that, since the exponent n 
is positive in the considered range, a lower i?„ = l/rf^ 
translates into an increased value of f2{R)' as r„ ap- 
proaches To ~ 12.8 Mpc from below, one eventually gets 
(R/Rn)" ^ 1 (in the inner, higher curvature region), in 
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FIG. 4. Plot of different dark matter mimicked density pro- 
files for different values of the characteristic length scale, 
r„ = 1.00, 6.31, 10.00, 15.84, 25.12 Mpc, where the dashing 
increases with the value of r„, compared with the observed 
dark matter component (black). These values are obtained 
for the best fit value n = 0.43. 



FIG. 6. Variation of 6(r„ = 10"^ pc, n) 




r(kpc) 



FIG. 5. Additional radial velocity resulting from the non- 
geodesic motion due to the dark matter component (dashed) 
and to the mimicked dark matter (dot-dashed) for different 
values of n, the same of Fig. |T|, and also for the total density 
profile (black). 



contradiction with the assumption of a perturbative non- 
minimal coupling. 



1. Non-geodesic motion 

In Fig. ([5]) one presents the plot of the additional ra- 
dial velocity Av — y/rf^ arising from non-geodesic mo- 
tion as a function of distance to the cluster's center, for 
the considered range of values for n. Notice that, de- 
spite the fact that the observable quantity is the velocity 
dispersion cr, estimated from the velocities of individual 
clusters, these curves enable one to see the relevance of 
the extra forces that appear due to the presence of the 
non-minimal coupling. 

Since the scalar curvature (i.e. the total density) is de- 




FIG. 7. Variation of fe(r„, n = 0.43). 



creasing, one concludes that the related force counteracts 
the gravitational attraction, as discussed in section p^II B| 
Thus, stability requires that this additional velocity be 
much smaller than the Newtonian radial velocity profile: 
this is attained for higher values of n. In particular, one 
concludes that a power-law coupling with n — 0.43 can 
effectively mimic the dark matter component of the Abell 
cluster A586, with a negligible effect on the velocity pro- 
file v(r). 

As already discussed in section |IIID[ one expects a 
decrease of the effect of the extra force for increasing ex- 
ponent n (for fixed r„). Indeed, as Fig. [s] depicts, the 
considered range always leads to a perturbative f2{R)- 
However, as Fig. |6] shows, values of n < 0.23 do not en- 
able a sufficiently strong suppression of the extra force: 
although the additional (outward) radial velocity is in- 
deed smaller than c, it may still be above the Newtonian 
counterpart. 

and the relevant values obtained 
one has e — 0.3626 and pdmo ~ 
ro « 12.8 Mpc. Inserting this into Eq. 

= 10"" 



Using Eq. 
from Ref. 
lO^^Mo/Mpc 



11) 



(41), together with Tc = 67 kpc and 



pc, one 

finds that the extra force is smaller than the Newtonian 
force when 
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&(10-2 pc,n) < 1 ^ n > 0.246 



(46) 



which should be compared with the aforementioned lower 
bound n > 0.23, obtained from Fig. [5] 

The lack of sensitivity of the mimicked dark matter 
density profile on the value of r„ may be used to establish 
a more stringent upper bound on this parameter. To do 
so, one fixes the best fit value n = 0.43 and searches for 
the value of r„ that fulfills the above condition: 



6(r„,0.43) < 1 



< 51 pc 



Variation of 6(r„ = 10 ^ pc,n) and 6(r„,n 
shown in Figs. [6] and [7] respectively. 



(47) 
0.43) are 



C. Relevance in other contexts 

The non-minimal coupling here considered may be re- 
garded as an approximation to a more evolved form for 
/2(i?), which would encompass couplings considered in 
previous studies, which perhaps hints at a Laurent se- 
ries, 



f2{R)= E 



R 

Rn 



(48) 



Thus, one is left with the possible conflict between as- 
trophysical relevant terms. If one considers three power- 
law couplings with n = —1, —1/3 and 0.43, the result- 
ing dark matter profile will not be a sum of the in- 
dividual contributions, due to the non-linearity of Eq. 
(23). Even if one disregards the resulting interaction be- 
tween terms (a reasonable approximation, as shown in 
Ref. [in]), the exponents simply yield the slope of each 
derived dark matter contribution, not their relative val- 
ues. Indeed, when numerically integrating Eq. (31), one 



imposes boundary conditions on the overall dark matter 
density {i.e. the dimensionless function g) — the rel- 
ative contributions of each dark matter component are 
then determined dynamically. 



Given the scaling relation Eq. (24), a larger n < 1 im- 
plies a steeper dark matter component Pdm{n) ~ 
Resorting to the p ^ r^^ outer Hernquist profile for the 
matter distribution of galaxies, one finds 



Pdmi-i) ~ r 



Pdm{0A3) 



Pdm(-l/3) 
„-7 



(51) 



Thus, the n — 0.43 merely gives rise to a quickly decaying 
component: the n = —1 and n = —1/3 contributions 
dominate at large distances, even if the former dominates 
at short ones. Conversely, the p ^ r~^^ behaviour of the 
cluster mass density translates into 



Each term would be valid in a particular regime: early 
vs. late time, central vs. long range, galactic vs. cluster 
size, etc. In particular, this series should also include the 
aforementioned n — —1 and n = —1/3 terms, relevant to 
mimic dark matter at galactic level, as well as a linear 
coupling R/Ri considered in the context of preheating 
and inflationary dynamics; the corresponding parameters 
are 



r-i 



{R-iY 

{R-l/3 



1/2 



21.5 Gpc 

1/2 _ -, «n „ -,n6 



(49) 



ri = i?i = 4.8VC X 10 



1.69 X 10" Gpc 

-29 



m 



with 10 < C < 10** [22]. 

Since the coupling terms relevant at galactic scales 
have negative exponents, it is trivial to conclude that 
their effect is negligible during inflation, when the curva- 
ture is much higher. Conversely, the effect of the terms 
involved in preheating and inflation are irrelevant at as- 
trophysical or cosmological scales at present. 

Similarly, the cosmological effect of the best fit cou- 
pling with n = 0.43 may be disregarded, as it is of the 
order of 



Rr 



R„ 



rn 



2n 



< 10 



-3 



(50) 



Pdm(-l) ^ r 



Pdm{0A3) r 



Pdm(-l/3) 
2.7 



(52) 



In order for the contribution from the n = — 1 and 
— 1/3 to be neglected, it suffices that they are subdomi- 
nant at short distances. One does not pursue an explicit 
computation of the effect of these three power-law terms 
here, given that the main scope was to describe the A586 
cluster. Nevertheless, it is not difficult to see that the 
difference between the typical galactic and cluster scales 
do allow for a consistent disentangling of the effects of 
these two sets of terms. 



V. CONCLUSIONS AND OUTLOOK 

In this work, the possibility of mimicking the dark 
matter component of a galaxy cluster in the context of 
a model, Eq. ( |13p , where the scalar curvature is non- 
minimally coupled with matter is discussed. For that, 
the cluster A586 has been examined, since it is relaxed 
and quality data for most of its relevant dynamical pa- 
rameters is available. 

For fi{R) = R, the best fit obtained with a power- 
law couphng f2{R) = (R/Rn)^ is for n — 0.43, with 
a strong constraint on the typical length scale r„ = 
l/\/Rn < 51 pc arising from the requirement that the 
induced extra-force is smaller than the Newtonian one. 
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A weaker constraint r„ ^ 



12.8 Mpc arises from the 



behaviour of the numerical integration of Eq. (31 ) 



The obtained fit does not correspond to the analytical 
value for n computed from the scaling relation between 
dark and visible matter, as the latter would give rise to 
a mimicked dark matter density that does not follow the 
observational curve (albeit it falls with the same slope). 
Instead, a shift from this analytical value is enforced, 
whereas the mimicked density is slightly steeper than the 
observed profile, but with a negligible difference in total 
mass. 

Furthermore, the exponent n is high enough so that the 
non- minimal coupling is perturbative, with f2{R) ~ 10^* 
suppressing the destabilizing effect of the outward extra- 
force arising from the non-conservation of the energy- 
momentum tensor, as previously found in the context 
of galaxy rotation curves (where two negative exponents 
n = — 1 and n = —1/3 were shown to correspond to the 



best fit for the rotation curves) [TU]. 

Studies addressing different scenarios have resorted to 
power-law couplings with different exponents than the 
best fit value n = 0.43 here reported — which may be 
regarded as another term in a putative series expansion 
of a more general non-minimal coupling. By the same 
token, yet unprobed terms of this series may be charted 
through the study of other phenomena and environments, 
where distinct curvatures and densities are at play. 

This first attempt at modelling the dark component 
of clusters shows that a consistent approach may be 
obtained by considering a suitable perturbative non- 
minimal coupling. In a subsequent work, this method will 
be extended to more dynamical clusters, using the avail- 
able profiles for the visible matter density and tempera- 
ture obtained from modifications of the usual /3-model, 
Eq. (|6]), and the polytropic temperature profile, Eq. ([7|, 
(cf. Ref. [33]). 
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